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ABSTRACT 

Aeroelastic stability analyses have been performed 
for the M0D-5A blade/aileron system. Various 
configurations having different aileron torsional 
stiffness, mass unbalance, and control system 
damping have been investigated. The analysis was 
conducted using a code recently developed by the 
General Electric Company - AILSTAB. The code 
extracts eigenvalues for a three degree of freedom 
system, consisting of: (1) a blade flapwise mode, 

(2) a blade torsional mode, and (3) an aileron 
torsional mode. Mode shapes are supplied as input 
and the aileron can be specified over an arbitrary 
length of the blade span. Quasi-steady aerodynamic 
strip theory is used to compute aerodynamic deriva- 
tives of the wing-aileron combination as a function 
of spanwise position. Equations of motion are 
summarized herein. The program provides rotating 
blade stability boundaries for torsional divergence, 
classical flutter (bending/torsion) and wing/aileron 
flutter. It has been checked out against fixed-wing 
results published by Theodorsen and Garrick. 

The M0D-5A system is stable with respect to diver- 
gence and classical flutter for all practical rotor 
speeds. Aileron torsional stiffness must exceed a 
minimum critical value to prevent aileron flutter. 
The nominal control system stiffness greatly exceeds 
this minimum during normal operation. The basic 
system, however, is unstable for the case of a free 
(or floating) aileron. The instability can be 
removed either by the addition of torsional damping 
or mass-balancing the ailerons. 

The M0D-5A design was performed by the General 
Electric Company, Advanced Energy Program Department 
under Contract DEN3-153 with NASA Lewis Research 
Center and sponsored by the Department of Energy. 
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INTRODUCTION 

Although aileron systems have widespread use on 
fixed-wing aircraft very few rotors have been 
designed with aileron controls. Large wind 
turbines, in particular, have used pitchable blade 
sections for power regulation and to start-up and 
shut-down. General Electric's 400 ft. diameter, 
7.3MW M0D-5A was originally designed with a pitch- 
able tip spanning the outer 25% of blade radius. As 
the design progressed, further studies determined 
that substantial weight and cost savings could be 
obtained by switching to an aileron control system. 
This provided the impetus for the work described in 
this work. 

In this paper we first describe the development of 
AILSTAB, a three degree of freedom stability 
analysis program. The results of the M0D-5A rotor 
blade stability analysis are then presented. Also 
included are the results of investigations, which 
are parametric in nature and show trends which 
should be similar for other WTG's. 

NOMENCLATURE 

a - distance, midchord to elastic axis, as 

percent of chord 

a - lift curve slope; lift coefficient per 

rad ian 

b - semichord 

c - distance, midchord to aileron hinge, as 

percent of chord 

c(k) - Theodorsen' s coefficient 

Ca - elemental aerodynamic damping matrix 

Ca - integrated aerodynamic damping matrix 

C$ - elemental structural damping matrix 
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r$ - integrated structural damping matrix 

e - distance, midchord to leading edge of 

aileron, as percent of chord 

h - blade flap bending degree of freedom 

la - torsional moment of Inertia of blade, less 

aileron, about elastic axis; per unit 
length 

Ig - torsional moment of inertia of aileron, 

about hinge; per unit length 

K h - elemental aerodynamic stiffness matrix 

1T A - integrated aerodynamic stiffness matrix 

Kcf - elemental centrifugal stiffness matrax 

1 Tqp - integrated centrifugal stiffness matrix 

Ks - elemental structural stiffness matrix 

K$ - integrated structural stiffness matrix 

Nb - blade mass, less aileron; per unit length 

Mg - aileron mass; per unit length 

M A - elemental aerodynamic mass matrix 


a 

rotor speed, radians/second 


0)f 

flutter frequency, radians/second 

U)h 

blade bending frequency, 
second 

radians/ 

(O a 

blade torsional frequency, 
second 

radians/ 

W B 

aileron torsional frequency 
second 

, radians/ 


METHODS OF ANALYSIS 



The AIlSTAB rotor blade stability analysis program was 
developed in a manner very similar to that which would 
be used for a fixed wing. The differences between 
rotor and fixed wing analyses are the variation of 
local velocity with span on a rotor blade, and the 
variation of stiffness with rotor rpm due to the 
centrifugal forces. The AILSTAB computer code can be 
used to predict divergence and classical blade bending/ 
torsion flutter, as well as aileron torsion/blade 
bending flutter. 
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integrated aerodynamic mass matrix 

elemental structural mass matrix 
integrated structural mass matrix 
blade radial station, dimensional 
Southwell coefficient, blade bending 
local velocity 

blade torsion degree of freedom 

aileron torsion degree of freedom 

critical damping ratio, blade bend- 
ing 

critical damping ratio, blade torsion 

critical damping ratio, aileron 
torsion 

air mass density; per unit length 

static moment of blade, less aileron, 
about elastic axis; per unit length 

static moment of aileron, about hinge; 
per unit length 

flapwise deflection mode shape 
blade torsion mode shape 
aileron torsion mode shape 
flapwise rotation mode shape 


SYSTEM DESCRIPTION 

The three degrees of freedom (DOF) in the analysis are 
blade flapwise bending (h), blade torsion (a), and 
aileron torsion (8). Figure 1 depicts these DOF and 
their sign conventions. The conventions are such that 
h is negative for a bending deflection toward the 
suction side of the airfoil, a is positive for a 
"nose up" rotation, and 8 is positive for an "aileron 
down" rotation and is measured relative to a. 

Three other parameters required for the analysis are 
depicted in this figure. All three are measured from 
the airfoil's midchord, are positive toward the trail- 
ing edge, and are expressed as a percentage of the 
semichord. The distance to the elastic axis is 
denoted "a" , the distance to the aileron leading edge 
is "e", and the distance to the aileron hinge Is 
denoted "c". 

Figure 2 shows the three DOF (mode shapes) depicted In 
three dimensions. 

ASSUMPTIONS 

The following set of assumptions, all of which are 
believed to be reasonable, were made in developing the 
computer code. 
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1. The equations of motion were linearized. 

2. Three degrees of freedom at a time, one 
flapwise, plus the blade and aileron torsion 
modes, are sufficient to determine the 
stability. 

3. Aerodynamic strip theory with no stall was 
used, i.e. the aerodynamic derivatives are 
independent of a- 

4. The local velocity is equal to fir, the 
rotational velocity times the radial 
distance, I.e. the free wind velocity is 
neglected. 

5. The Theodorsen coefficient, c(k), is equal 
to 1.0, i.e. Quasi-steady aerodynamics are 
used. This should give conservative results 
for both blade-bending/torsion flutter and 
blade-bending/aileron torsion flutter. 

6. Aerodynamic derivatives for an unsealed gap 
(ref. 1) are used if c£e. 


EQUATIONS OF MOTION 

The equations of motion were developed for a 
representative airfoil element of length “dr" and 
Integrated along the span of blade with weighting as 
determined by the mode shapes. The aerodynamic 
equations incorporated in this analysis were those 
of Smilg and Wasserman (ref. 1). Inertial equations 
of motion were derived with centrifugal stiffening 
terms added. The final form of the equations is: 


D<S*KcF- k a] ' 


where the matrices subscripted S (structural ) and A 
(aerodynamic) are composed of elemental mass, 
damping, and stiffnes terms integrated along the 
blade span with modal weighting. 


py -J L0 h ** 0 0 ] [M S ] 
span 



A detailed description of the terms in the elemental 
matrices is presented in the appendix. 

SYSTEM STABILITY 

In order to determine the blade's flutter stability, 
the integrated mass, stiffness and damping matrices 
are formed into a six by six dynamical matrix from 
which complex eigenvalues and eigenvectors are 
determined. The form of the dynamical matrix is: 



where 

M 3 -M A + 

C 3 -C A + C S 

K 3 -K A + Ks + KcF 

I *3x3 unit matrix 

The critical damping ratios (c), and the 

frequencies in hz (f^)» are determined from the 
eigenvalue (R) as follows. 

5 = - REAL(R) , f f - ABS(R) 

ABS(R) 2 7T 

The output of the AILSTAB stability analysis program 
is eigenvalues, and eigenvectors if desired. The 
program is organized so that a series of cases may 
be run for a particular configuration, with rpm 
varied. The critical damping ratio and coupled 
frequencies are determined from the complex eigen- 
values, and the damping in each mode can be plotted 
vs. rpm to Illustrate system stability. In the case 
of blade bending/aileron tension flutter there is a 
range of rpm between which the instability exists. 
By plotting the range of unstable rpm vs. a design 
parameter such as aileron control system stiffness, 
aileron damper, or mass balance, a stability 
boundary may be constructed. 


The stiffness contribution due to centrifugal 
stiffening, K^p, is formed similarly with a mode 
shape of flapwise rotations substituted for the 
flapwise deflection mode shape. 



DESCRIPTION OF ANALYZED SYSTEM 

The M0D-5A is a 7.2MW wind turbine with a teetered 
rotor. Ailerons on the outer 40% of the 200 ft. 
radius blades are used to regulate power and to shut 
down. The ailerons are hinged at their leading edge 
and are 40% of the chord width. 

Three blade flapwise mode shapes were used in the 
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analysis- They were 1) the teeter mode with a 
frequency of 1 per rev, 2) the 1st collective with a 
frequency of 7 radians/second, and 3) the 1st cyclic 
with a frequency of 13.8 radians/second. These mode 
shapes were calculated for an isolated blade ( 1 . e - , 
not attached to the wind turbine). The collective 
mode of the isolated blade is found by providing a 
cantilevered root condition In the flapwise 
direction. The cyclic mode is determined by 
providing a pinned root condition in the flapwise 
direction. A plot of these three flapwise modes Is 
presented in Figure 3. Southwell coefficients may 
be input to the program so that both the collective 
and cyclic frequencies may be varied with rpm to 
account for the varying centrifugal stiffening. For 
the M0D-5A analysis the important instability 
occurred at a low enough rpm so that the centrifugal 
stiffening was not important to the results. 

The three aforementioned flapwise modes were each in 
turn analyzed in combination with the blade torsion 
mode shape and an aileron torsion mode. Higher 

modes than these were also analyzed, but were not 
found to be critical. The blade torsion mode had a 
frequency of 51 radians/second. By comparison, the 
ailerons are essentially rigid in torsion with 
cantilevered frequencies above 400 radians/second. 
For all practical purposes, the aileron natural 
frequencies are dominated by the control system 
stiffness and oscillate as a rigid body. Rather 
than attempting to model the actual aileron tor- 
sional natural mode, the frequency, or equivalently 
the actuator stiffness was varied, to determine the 
minimum requirements. In this way failure modes, 
such as loss of actuator hydraulic stiffness, are 
fall-outs of the analysis. In addition to aileron 
frequency sweeps, variations in aileron torsional 
damping, mass-balancing, and aileron spanwise length 
were considered. The ailerons center of gravity is 
aft of the 60% chord hinge line. The aft center of 
gravity has a de-stabil izing effect. 

ANALYSIS RESULTS 

The most critical condition will be discussed first. 
It occurs when the root torsional stiffness provided 
by the actuators is lost and the aileron is free to 
rotate about its hinges. This cannot happen under 


normal circumstances, so it represents a system 
failure. Stability boundaries are presented in 
terms of control system stiffness, control system 
damping, and the degree of mass-balance. A final 
case considers the stability of an aileron spanning 
only the outer 27.5% of the blade. 

FREE AILERON 

Figures 4a-c show damping vs. rotor speed for the 
baseline blade with a free, unbalanced aileron (the 
aileron torsion frequency of 1 per rev or IP, is due 
to centrifugal stiffening). Below each damping 
curve, the natural frequencies are plotted vs. rpm. 
Both coupled (dashed lines) and uncoupled (solid 
lines) frequencies are shown. At rpm's where 
uncoupled frequencies coincide, a decrease in sta- 
bility is noted in the corresponding damping curve. 
The Figures 4a, b, and c, illustrate the stability 
with the teeter, flap collective, and flap cyclic 
modes, respectively. Aileron torsion coupling is 
seen to cause an instability only with the flap 
collective mode. The ailerons are unstable in the 
region of low rotor speed, 3-12 rpm, which is 
typical of wing and aileron systems with an unbal- 
anced mass. In particular, there Is the possibility 
of instability when the aileron torsional frequency 
is less than the flapwise frequency. The system in 
Figure 4 becomes stable again at 12 rpm, because the 
torsional aerodynamic spring increases the aileron 
frequency beyond that of the first flapwise mode. 
The instabilities, which are seen in all plots 
between 55 and 60 rpm, are classical bending-torsion 
flutter of the blade. 

STABILITY BOUNDARIES 

A flutter boundary for the M00-5A blade with 
unbalanced aileron is given in Figure 5. To 
generate the boundary, the aileron root torsion 
spring was Increased in increments to find the 
stiffness at which the torsion mode became stable. 
At any value of stiffness where an instability 
occurred, the values of rpm between which the mode 
was unstable were found and plotted. This figure 
shows that an aileron torsional frequency of 7.5 
radians/second is needed to provide neutral stabil- 
ity. This same procedure was followed for the 
addition of aileron torsional damping rather than a 
spring. The resulting flutter boundary is presented 
in Figure 6. 
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The comparison of stiffness and damping requirements 
Is an interesting sidelight to the stability problem. 
If the damping rate is multiplied by the flutter 
frequency, the effective impedance, in stiffness 
units of the damper is found. Figure 7 contains 
plots of impedance vs. flutter damping ratio at 
6 rpm for both spring and damper systems. The 
system's stability Is largely a function of the 
aileron torsional impedance whether it be derived 
from a spring or a damper. This conclusion is 
further strengthened by Figure 8 which shows the 
stability boundaries in terms of impedances. The 
approximate equivalence of spring and damper 
impedance effects is an important consideration 
during dynamic conditions, such as pitch change in 
which the hydraulic actuator impedance has both 
spring and damper characteristics. 

To prevent flutter without need for a mimimum 
aileron torsional stiffness or damper, balance 
weights would have to be added to the ailerons. 
With the ailerons unbalanced, the minimum damping 
ratio calculated in the AILSTAB rpm sweep was 
approximately -12%, as can be discerned from Figure 
4b, The variation of modal damping with RPM is 
shown for a fully (100%) mass-balanced aileron in 
Figure 9. The system is stable. The variation of 
minimum damping in the aileron mode is shown for 
varying degrees of mass-balance in Figure 10. 
Neutral stability can be obtained with an 85% 
mass-balanced system. 

EFFECT OF AILERON LENGTH 

Similar analyses to those discussed above were 
performed with the free aileron section extending 
from .725 radius to the tip, rather than from .60 
radius. In this configuration the different modal 
weighting caused an instability of the aileron 
coupling with the blade cyclic bending mode. Aileron 
torsion coupling with the blade collective bending 
mode also produced an instability, as it had with 
the longer aileron. 

Since the shorter aileron was unstable in coupling 
with the higher frequency cyclic flapwise mode with 
a flutter frequency of approximately 14 radians/ 
second, a higher dimensional damping coefficient was 
required to stabilize it. The longer aileron had 
unstable coupling only with the collective mode, 


which had a flutter frequency of approximately 7 
radians/second. 

The damping vs. rpm plots for the 27.5% span, free, 
unbalanced aileron analysis are presented in Figure 
11 . 

Stability boundaries of rpm vs. aileron frequency 
are plotted in Figure 12 to show the effect of added 
root torsional stiffness. Boundaries for the addi- 
tion of torsional damping are shown in Figure 13. 
To again demonstrate the similarity of results from 
adding impedance, whether from stiffness or damping, 
rpm vs. impedance stability boundaries are shown in 
Figure 14. 

CONCLUDING REMARKS 

The free unbalanced aileron caused the system to 
become unstable either with a length of 40% or of 
27.5% of blade radius. 

These instabilities can be removed with the addition 
of Impedance to the aileron torsion degree of 
freedom. The actuator stiffness normally supplies 
an impedance well in excess of that required, but on 
the M0D-5A torsional dampers have been added to 
protect the system in the event of an actuator 
system failure. These dampers are passive elements 
which will always be operative. The damper forces 
far enough below those which are present due to the 
aerodynamic forces in normal operation so that their 
presence will not penalize control system design. 

An alternate method of stabilizing the system would 
be through the addition of balance weights to the 
aileron. This method was deemed unwieldy and tor- 
sional dampers chosen instead. 
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APPENDIX 

Equations of Motion 


southwell coefficient for the flapwlse mode. It is 
an approximation used to relate the rotating and 
non-rotating blade natural frequencies. 




to 2 

ROT 

■ “ 2 NON- 

k cf 

(1.2) 

= 

r a 2 

: M s (1,2) 

k cf 

(1.3) 

= 

r a 2 

1 M s (1 ,3) 

k cf 

(2.3) 

* 

fl z Me (2,3) 

k cf 

(2,1) 


k cf 

(1,2) 

k cf 

(3,1) 

= 

k cf 

(1,3) 

k cf 

(3,2) 

= 

k cf 

(2,3) 


STRUCTURAL (Left Hand Side) 



+ [C S ] 


r 


L S J 


+ ^ 



(1,1) » M a + Mjj 

^(1,2) • a a + a B + b ^ C ‘ a ^ M B 

Ms(l,3) » a 3 

M s (2,2) « I a + I fi + 2b (c-a) 0 g + b 2 (c-a)* Mg 

M s (2,3) = I 0 + Ig + b (c-a) ag 

M s (3,3) = Ig 


Other structural mass terms are symmetric. 


CjO.l) - 

Cj (2,2) ■ 

C S t 3,3 ) = ^6 ? 0 

Other damping terms are zero. 

K s (1.1) * o> h 2 M a + S h fl* M s (1,1) 

K S < 2 > 2 ) * “a" J a + M S ( 2 > 2 ) 

K$ (3,3) = cog 2 I B + fl 2 M s (3,3) 

Other stiffness terms are zero. 

The above structural mass, stiffness, and damping 
matrices are all multiplied by mode shapes at each 
radial station and integrated. 

i .e. 

N T 

M S = l [0h 0 O 00] [M S ] [0h 003 Ar 

0 

S h in the above stiffness equation, is the 


Diagonal terms are zero. 


This centrifugal stiffness matrix is multiplied by 
mode shapes at each radial station and integrated. 
Unlike the structural stiffness matrix, the first 
row and column are multiplied by the modal rotation 
rather than deflection. 

i .e 

N T 

k CF ■ l [00h 9 a Ob] [K C f] [0Qh Qa 
0 

Where, (J 0h = d0 h 
dr 


AERODYNAMIC (Right Hand Side) 




[M A ] 




\ + C c a 3 i 


[k a ] 


m a 

0.1) 

- -7i p b 2 


m a 

(1,2) 

“ ir p ab 3 


m a 

0,3) 

' P (T, ♦ 

(c-e) 0 3 ) b 3 

m a 

(2,2) 

= -Ti p (a 2 

+ 1/8) b" 

m a 

(2,3) 

■ p [T? + 

(e-a) Tj + 1/4 (c-e) 



- (a h 

■ 1/2) (c-e) 0,] b' 

* 

(3,3) 

- t [T* 

+ (c-e) 037 - (c-e) 2 


ir 


Other mass terms are syimetric. 


C A (1.1) = -p a a r b c (k) 
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C A ( 1 . 2 ) “ 
C A ( 1 . 3 ) - 

C A (2,1 ) * 
C A (2,2) * 
C A (2,3) * 

C A (3,1) = 
C A (3,2) = 

C A (3,3) = 

K A (1,1) ■ 
K A (1,2) = 

Ka (1,3)) 

( 2 , 1 ) - 
K A (2,2) 

K A (2,3) 

K a (3,1) 

Ka (3,2) 

Ka (3,3) 


-1/2 paarb 2 +pa(a- 1/2) a r b 2 c (k) 

1/2 P_ ¥ T. n r b 2 + P a [ (c-e) 0i - 1/2 Tn] a r b 2 c (k) 

IT IT 

p a (a + 1/2) a r b 2 c (k) 

1/2 p a (a - 1/2) a r b 3 + p a (1/4 - a 2 ) ft r b» c (k) 

P ¥ [P - 1/2 (a - 1/2) T*. + 1/2 (c-e) 0>] arb’ 

TT 


+ P a [(a + 1/2) Tn - (a + 1/2) (c-e) 0i] arb’c (k) 
n 2 


P a [(c-e) - 1/2 Tm] arb'c (k) 


TT 

1/2 P a [(c-e) 032 - (P - Ti - 1/2 T-)] a r b 3 

TT 

+ P a [1/2 (a - 1/2) Tm + (1/2 - a) (c-e) 03i] a r b 3 c (k) 

TT 

1/2 P a [1/2 T* Tn + (c-e) (03« + 0m) - (c-e) 2 0a»] a r b 3 
i* 

+ 1/2 P "a [-1/2 Tn Tm + (c-e) (02 03i + 0i 0.) 


-2 (c-e) 2 0i 03i] a r b 3 c (k) 


0 

-p a" (a r) 2 b c (k) 

= - £ ¥ Tio (si r) 2 b c (k) 

TT 


p a (a + 1/2) (Q r) 2 b 2 c (k) 

-1/2 P_ ¥ [T. + Tm] (a r) 2 b 2 + £ ¥ (a + 1/2) Tio (a r) 2 b 2 c (k) 
it it 

0 

+ P ¥ ((c-e) 03i - 1/2 Ti 2 ) (a r) 2 b 2 c (k) 

IT 

1/2 p a [(c-e) 035 - (Ts - Tt Tio)] (a r) 2 b 2 
¥ 2 

+ 1/2 P 7 [2 (c-e) 0i 03i - Tio Tm] (a r) 2 b 2 c (k) 
P 
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Ti 


T 3 

T* 

T s 

h 

Ti o 
hi 

T n 
P 

0» 

0 2 

03 

05 

06 
08 
03 1 
01 0 
017 
03 2 
03 5 

03 6 
03 7 


= -1/3 sin (cos" l C) (2+C 2 ) + C COS-'C 
= -(1/8+C 2 ) ( COS" 1 C ) z + 1/4 C SIN (COS - 1 C } COS-'C (7+2C 2 ) 

* -COS-'C + C SIN (COS-‘C) 

* -1 + C 2 - (COS - 1 C ) 2 + 2C SIN (COS- 2 C) COS-*C 

■ -(1/8+C 2 ) C0S-‘C + 1/8 C SIN ( COS - 1 C ) (7+2C 2 ) 

* SIN (COS-^C) + COS-‘C 

* (COS" 1 C ) C ( 1-2C) + SIN (COS'*C) (2-C) 

« SIN (COS" 1 C) (2+C) - COS'«C (2C+1) 

= -( SIN (COS' 1 C) ) 2 1/3 

= n - COS' 1 (-e) + SIN ( CO$- l (-e) ) 

« ( ir-COS'U-e) ) ( l-2e) + SIN ( COS'M-e) ) (2-e) 

= ti - COS'M-e) - SIN ( COS'M-e) ) e 
= SIN ( COS'M-e) ) (1+e) 

= 2 ( Tr-COS'M-e) ) + SIN ( COS'M-e) ) 2/3 (2+e) (l-2e) 

= ( ir-COS'M-e) ) (-l-2e) + SIN ( COS'M-e) ) (2+e) 

= ir -COS'M-e) - SIN ( COS'* (-e) ) 

1 03 1 05 

■= 0 3 2 + [SIN (COS'M-e)]" 

= TT -COS'M-e) + SIN ( COS'M-e) ) (l+2e) 

* 2 [SIN ( COS'M-e) )] 2 

= 032 03 + 2 [SIN (COS'M-e) )]" 

= 03 (02 - 03 ) 


The aerodynamic mass, damping, and stiffness matrices 
are all multiplied by their mode shapes at each 
radial station and integrated. 

i .e. 

- N T 

m A = l [0p 0 a 0 b1 C m aJ [0h 0 a 00] Ar 

0 

DIVERGENCE 

Torsional divergence, if present, will show up in 

the roots of the stability equations. The following 
has been added so that the divergence speed, which 
often lies beyond the RPM range of interest, may be 
computed directly. 

To determine the blade's divergence speed the square 
of the rotational rotor velocity, q 2 , must be 
factored out of the lower-right-hand 2x2 partition 
of the integrated aerodynamic stiffness matrix, 


K A * The same partition is factored out of the 
structural stiffness matrix i< s (the Southwell 
coefficient terms were ignored here). The two 
resulting partitions are then set equal to each 
other and the characteristic equation solved as 
shown below. 


K A ' = 

*a 

(2,2) 

Ka 

(2,3) 

1 


jj<A 

(3,2) 

Ka 

(3,3) 

fl 2 

«S' ' | 


(2,2) 


(2,3) 



3 

(3,2) 

% 

(3,3) 


KS‘(1,1) - 

Ka‘ 

(i,D a 

«s 

MM2) 

- K A '(1,2) ft 2 

Ks' (2,1) - 

•V 

(2,1) ft 2 

KS 

M2, 2) 

- Ka' (2.2) a 2 

K S M2,1) = 

K 

S M1,2) 

= 0 
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Expanded, the resulting equation is: 


a ( k A * ( 1 » 1 ) K A * (2,2) - K a '(1,2) K a '(2,1) ) 

(K s , (2,2) K A '(1,1) + K s '(l,l) K A '(2,2) ) + K' s (l,l) K' s (2,2) = 0. 

which is easily solved fori), the flutter speed in radians/second. 



a = DISTANCE - MIDCHORD TO ELASTIC AXIS 
e = DISTANCE - MIDCHORD TO AILERON L/E 
c = DISTANCE - MIDCHORD TO AILERON HINGE 
h = BLADE FLAPWISE DEFLECTION 
a = BLADE TORSION ANGLE 
6 * AILERON TORSION ANGLE w.r.t. 


M'y 



L' = RUNNING LIFT FORCE 

- POSITIVE POINTING FROM SUCTION 
SIDE TOWARD PRESSURE SIDE AS DEFINED 
IN REFERENCE 1 

M'y = RUNNING MOMENT ABOUT ELASTIC AXIS 
M' h = RUNNING AILERON HINGE MOMENT 


Figure 1. Sign Conventions and Terminology 
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Figure 4c. Flap Cyclic Mode Stability 
Free Unbalanced Aileron 



AILERON FREQUENCY 
(rad/sec) 


DAMPING COEFFICIENT 
(Ib-ft-sec) 


r !gure 5. MOD- 5 A Flutter Boundary - Shown in 
Terms of Aileron Torsional Frequency 


Figure 6. M0D-5A Flutter Boundary 
Terms of Aileron Hinge 


r~ 

4000 


- Shown in 
Damping 


no 



Variation in ftodal Damping with Aileron 
Torsional Impedance, Shown for 6 RPM. 


stability of Free 100% Balanced Aileron 
60 - 1 .OR 


CRITICAL DAMPING RATIO CRITICAL 0AMP1NG RATIO CRITICAL DAMPING RATIO 
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COLLECTIVE MODE 
fn * 7 rad/sec 


RPW 


STABLE 

REGION 


10 




AILERON TORSION FREQUENCY 
(rad/sec) 


AILERON TORSION FREOUENCY 
(rad/sec) 


Figure 12. N0D-5A Flutter Boundaries 

Unbalanced Aileron . 725 R - 1 .OR 
Shown as a Function of Aileron 
Torsional Frequency 


COLLECTIVE NODE 
fn * 7 r/s 



CYCLIC NODE 
fn = 13.8 r/s 



Figure 13. f^OD- 5A Flutter Boundaries 

Unbalanced Aileron _ 725 R - 1 .OP. 
Shown as a Function of Aileron 
Damping 
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COLLECTIVE MODE 
fn * 7 rad/sec 


CYCLIC MODE 
fn * 13.8 rad/sec 




AILERON TORSIONAL IMPEDANCE AILERON TORSIONAL IMPEDANCE 

10 3 ft-lb/rad 10 3 ft-lb/rad 


Figure 14. MO P-5 A Flutter Boundaries 

Unbalanced Ailerons .725R - 1 ,0P. 
Shown as a Function of Aileron 
Torsional Impedance 
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